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We discuss the implications of recent experimental evidence of nuclear Short Range Correlations 
(SRCs) on the modeling of the wave function of complex nuclei. We perform a calculation of 
potential energy contributions of pp and pn pairs in nuclei, showing that the presence of strong 
tensor correlations produce a ratio between the twos of about 1/9, which strongly deviates from the 
^\ . combinatorial counting of the number of pairs. We also discuss implications for the production of 

' Monte Carlo configuration for the simulation of nuclear reactions involving complex nuclei. 

o : 

. In the past decade, various measurements have been unambigously identified two-nucleon Short Range Correlations, 
Q ' studied their structure and related them to the underlying basic short range Nucleon-Nucleon (NN) interaction [H-[l]- 
/~N \ In [l| two- and three-body correlations were studied in large nuclei and related to deuteron and ^He measurements; 
' the (p, 2p-\-n) 0] experiment studied the directional correlation between proton and neutron momenta, while in f3| it 
was shown that about 90% of high momentum protons are correlated with a neutron; in ^2|| high momentum protons 
• in (e, e'p), (e, e'pp) and (e, e'pn) were investigated, finding a recoling, back-to-back proton in 10% of the events and a 
neutron in 90% of the events, consistently with After these experiment, several nuclear theory groups showed that 
, the measured ratio is a strong indication of the operation of an NN tensor force in the pair at the nucleon separations 
I ' and relative momenta studied d, d. Recently, an effective many-body approach for the description of ground-state 
'■^ , wave functions have been adopted for complex nuclei by the Perugia group, relying on a cluster expansion for the 
^ ' calculation of expectation values of quantum operators on the ground state wave functions with realistic forces 0] . In 
(~| this work, we investigated the effect of correlation on the fraction of potential energy in the nucleus which is carried by 
the different pn and nn pairs; moreover, we investigated the possibility of inclusion of full state-dependent correlations 
in nuclear configurations for the simulation of NA and AA collision at high energies. 
" ^ ' We can calculate the potential contribution to the ground state energy according to 



O : (^) = E(*(^^^)) = E(E-^"^(-^^)0?^) = ^^^^ E / dr,dr,p(:i\r,,r,)v^-\r,,). (1) 

(2) 

Here pn ^2) is the state-dependent two-body density matrix 

r ^ 

P^^\rur2) = / \{dr,r{ri,.:,rA)d^;$i^{r^,...,rA), (2) 



which we have evaluated within the cluster expansion method with ?A(ri, ...,rA) — Y[t<j /(''y O*/*!^!; being 

the independent particle model wave function, and is the operator 

Ou e {i,(Ti .^2,^12} {i,xi-T2} , (3) 

acting between particles 1 and 2, present in both the nucleon- nucleon potential and in the ground state wave function. 
The two-body density appearing in Eq. ([1]) can be splitted within our approach into the contributions due to proton- 
proton, proton-neutron and neutron-neutron pairs and the potential energy can thus be splitted into the corresponding 
contributions. The radial two-body density 

P^n\ri2) = / dRp^^Uri^R+^ri2,r2^R-^rr2] (4) 



Hri2) = J dRp^^^ (^ri^R+ iri2 ,r2^R- ^rr2^ 



*Also at: Department of Physics, University of Perugia and INFN, Sezione di Perugia, Via A. Pascoli, 1-06123, Italy 



2 




is shown in Fig. [T] for ^^O and "^"Ca, for n = c,a,T,aT,S,ST in Eqs. ([2]), Q and ([4]). We have calculated the 
pp and pn contributions to potential energy for ^^O and '"'Ca, using the density of Eq. ^ and the Argonne AV8' 
potential, obtaining for the following results. The contribution from central correlations is small and in this case 
the probabilities from pp and pn pairs are exactly proportional to the fractions estimated combinatorically, namely 
Z{Z-1)/A{A-1) and 2ZN/A{A-1), respectively, giving P{pp) = 23% and P{pn) = 53% for ^^O. In the case of full 
correlation, this proportionality docs not hold anymore, and we find P{pp) — 8% of the total and P{pn) — 83% of the 
total, or P{I = 0) = 74% and P{I = 1) = 26%, / being the total isospin of the pair. The corresponding calculations 
for "^^Ca give practically the same results. 

When simulating the collision of two heavy ions, one starts by constructing some picture of the two involved nuclei, 
with the positions of the nucleons are sampled according to a density distribution function; this approach completely 
ignores the structure of the wave function of the nucleus, which is a highly complicated object depending on the 
positions, momenta, spin and isospin states of A nucleons. In Ref. [1] it was developed a Monte Carlo code to provide 
a more realistic implementation of the nuclear wave function, including central correlations, and providing a good 
description of two-body densities of the nucleus as compared with the ones obtained within an independent particle 
model. From this work, we know that basic quantities as the potential energy of the nucleus are strongly affected by 
realistic SRCs, therefore it appears mandatory to include state- dependent correlations in the approach of Ref. Q. To 
this end, the first improvement to the mentioned approach was to distinguish protons and nucleons, and to implement 
state-dependent correlations between first-neighbor nucleons. With the newly developed configurations 9] we will be 
able to investigate various aspects of AA collision Fig. [5] shows the spectator systems after a Pb — Pb collision, 
the NN scattering being treated within the Glauber theory. Correlated nucleons, which will be emitted with high 
momentum, are depicted in a different color. 
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FIG. 2: The spectator systems after a Pb — Ph collision along the z direction. Left: side view; right: view from behind. Red and 
blue spheres are protons and neutron, respectively, while the green ones are those nucleons which belonged to a correlated pair, 
before the interaction, their correlated partner being among the interacting nucleons (hidden, in this figure). The axes units 
are in fm and the dimension of the spheres are taken as the rms charge radius of the proton; Glauber parameters correspond 
to RHIC energies. Animations are available 9] along with the configurations used for the colliding nuclei. 
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